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ABSTRACT: We propose to split the sparticle spectrum from the hierarchy between the GUT 
scale and the Planck scale. A split supersymmetric model, which gives non-universal gaugino 
masses, is built with proper high dimensional operators in the framework of SO(10) GUT. 
Based on a calculation of two-loop beta functions for gauge couplings (taking into account all 
weak scale threshold corrections), we check the gauge coupling unification and dark matter 
constraints (relic density and direct detections). We find that our scenario can achieve the 
gauge coupling unification and satisfy the dark matter constraints in some part of parameter 
space. We also examine the sensitivity of the future XENONIT experiment and find that 
the currently allowed parameter space in our scenario can be covered for a neutralino dark 
matter below about 1.0 TeV. 
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1 Introduction 


The discovery of a 125 GeV Higgs boson [1, 2] by both the ATLAS and CMS collaborations 
has completed the Standard Model (SM). The experimental data reported so far by LHC 
agree quite well with the SM predictions. On the other hand, the lack of significant hints in 
electroweak precision tests and the absence of evidences for new particle contents at the LHC 
challenge many proposals of new physics beyond the SM, including weak scale supersymmetry 
(SUSY). 

SUSY, which was regarded for a long time as one of the most appealing extensions of 
the SM, has many desirable features. For example, the observed 125 GeV Higgs boson falls 
within the narrow 115—135 GeV window predicted by the minimal supersymmetric standard 
model (MSSM). Besides, the genuine unification of gauge couplings [3-6], which can hardly 
be achieved in the SM, can be successfully realized in the framework of low energy SUSY. 
Also, the puzzle of cosmic dark matter can naturally be explained in SUSY. 

Although SUSY is appealing, null search results of sparticles at LHC suggest that either 
low energy SUSY needs to be tuned or sparticles are well above the weak scale. In fact, the 
LHC data has already set a limit [7, 8] mg > 1.5 TeV for mg ~ mg and mg Z 1 TeV for 
mg >> mg for certain popular CMSSM models. On the other hand, the observed mass of 
the Higgs boson requires rather large loop effects of top squarks in the MSSM and CMSSM, 
which implies some extent of fine-tuning [9]. So naturalness in SUSY may be realized in a 
more involved way even though it was initially proposed to solve the hierarchy problem. 

Split SUSY, proposed in [10-12], gives up the naturalness criterion while keeps the other 
two main advantages: the gauge coupling unification and viable dark matter candidates. This 


scenario assumes a very high scalar mass scale Mg and the low energy spectrum contains only 
the gauginos and higgsinos as well as a fine-tuned Higgs boson in addition to the SM sector. 
The SUSY CP and flavor problems can naturally be solved in this scenario due to very heavy 
sfermions. The latest results of the Higgs mass from the LHC measurement, together with the 
requirement of gauge coupling unification, suggest a scalar superpartner mass scale roughly 
of order Ms ~ 100 — 1000 TeV [13-15], which indicates a moderately split spectrum for 
split SUSY. The hierarchy of sparticle spectrum is argued to be determined by the gauge 
loop factor [14]. We propose in this paper to split the sparticle spectrum from the hierarchy 
between the GUT scale and the Planck scale. 

Actually, there are various ways to split the sparticle spectrum. In this work we propose 
to use the generalized gravity mediation [16-18] with non-renormalizable Kahler potential 
and superpotential. Because GUT is one of the retained motivations of split SUSY, cer- 
tain high-representation Higgs fields of GUT group could appear in the non-renormalizable 
Kahler potential and superpotential. Such high-representation fields could not only amend 
the gauge coupling unification condition at the GUT scale, but also establish new relations 
among various theory inputs. In our scenario we have non-universal gaugino masses at the 
GUT scale (for other scenarios giving non-universal gaugino masses, see [19-21]). Such non- 
universal gaugino masses can naturally appear with non-renormalizable Kahler potential and 
superpotential involving various high-representation Higgs fields. 

This paper is organized as follows. In Sec. 2 we present our model. In Sec. 3 we 
check the phenomenology of our model. Based on a calculation of two-loop beta functions for 
gauge couplings (taking into account all weak scale threshold corrections), we check the gauge 
coupling unification and the dark matter constraints. The future XENONIT sensitivity to 
our scenario is also examined. Finally, in Sec. 4 we give our conclusions. 


2 A split SUSY model from SUSY GUT 


There are many possible ways to mediate the SUSY breaking effects from the hidden sector 
to the visible sector. A very interesting and predictive possibility is the gravity mediation. 
With certain non-renormalizable terms, proper soft SUSY breaking parameters can be gen- 
erated. In many popular gravity mediation scenarios, the Kahler potential is assumed to 
be minimal. However, a general Kahler potential seems to be more natural. When certain 
high-representation chiral fields for the GUT group are involved in the non-renormalizable 
Kahler potential, the kinetic terms of superfields could have other contributions after the 
GUT symmetry breaking. New non-renormalizable terms in the superpotential involving 
high-representation fields could also be important. 

In general, the non-vanishing F-term VEVs of certain fields which break SUSY could 
be either gauge singlets or non-singlets. In this section, we propose that SUSY breaking is 
triggered by the GUT group non-singlet F-term VEVs. The gaugino and sfermion masses 
will be generated by some non-renormalizable operators which could arise from integrating 
out certain gravitational effects. In order to get more compact spectrum and simplify the 


relevant expressions, we adopt the SO(10) GUT group with the Georgi-Glashw SU(5) GUT 
group as an intermediate stage in the symmetry breaking chain: 


SO(10) SU(5) x U(1)x SU (3). x SU (2), x U(1)y. 


16,16 54(45) 


Here U(1)x is broken by possible v? component VEVs of H16 and His Higgs fields. Therefore, 
most of the results in this paper are also valid in SU(5) GUT. The possibility of SO(10) to 
Pati-Salam will be considered elsewhere. 

In order to accommodate the non-minimal Kähler potential, we need to know the group 
products of various SO(10) representations [22]. The spinor representation (for the matter 
part) can be decomposed as 


I6 & 16 = 1 $ 45 ẹ 210, (2.1) 

while the adjoint representation for gaugino is 
(45 Q 45)symmetric = 1 © 54 $ 210 @ 770. (2.2) 

and the fundamental representation of Higgs is 
10 & 10 = 1 $ 45 $ 54. (2.3) 


We assume the non-minimal gauge kinetic term for vector supermultiplets 
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with M, being the reduced Planck scale. So the kinetic part for gauge field with non-minimal 
Kahler potential is given by 


1 n 
L= -gir EE] - TEM Tr|Fu bF], (2.5) 
with k being the normalization factor for various representations according to Tr(T?T°) = 
ko”. After the GUT non-singlet develops a VEV, (©) = v + Fe6? with v? > Fe, the 
unification condition turns into 
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with 631 being the appropriate group factors for SU (3)e, SU (2)z, U (L)y , respectively. If the 
hierarchy between the (first step Georgi-Glashow SU(5)) GUT scale and the Planck scale is 
not small, the previous GUT conditions turn into two independent new GUT conditions: 
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For non-singlet ©, the F-term Fg can be decomposed as (F6)ap = Fy - Aah with Aap being 
the group factor and Fy the universal part. The review of the group structure can be found 
n [22]. 

The group structure of the 24 component F-term VEV of 54 representation Higgs can 
be written in terms of 10 x 10 matrix 


3,1 1 1 1 1 1 1 1 1 1 
F; = FPA TE =, =, Sooo oa Oo) 2.9 
( 54) ab 54/4ab 54 5l 3’ 3’ 3’ 2’ 2? 3’ 3? 3° 9’ 5) ( ) 
The gaugino will get contribution from 
< Fọ > 
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with the F-term VEVs (F54)ab = FY, Aab. Here the universal part F, EA is independent of the 
group structure. Then the non-universal gaugino masses are given by 


= a Mə = A tis M3 = TET (2.11) 


FU 
my/2 = va (2.12) 


with 


The sfermion masses and kinetic term will be generated by the following non-renormalizable 
Kähler potential 


1(®'®) ado (2.13) 


K = -7p 


with proper F-term and lowest component VEVs of ®, respectively. Here we assume that 
the universal part of the kinetic terms ®'® is approximately canceled by a similar high 
dimensional operators with a lowest component VEV of singlet. 

We know from the group theory that possible contributions to Kahler potential for matter 
content (filled in 16 representation of SO(10)) can arise from the following type of Higgs fields 


1945554854, 16210516816. (2.14) 


So the F-term VEV of the 54 representation chiral superfield will contribute to non-universal 
sfermion masses of order (F54/M,)? in addition to possible contributions from Fie. We 
assume F16 ~ F54 and we always have Fyg = Fig. Since vie > v54, the 54 representation 
field will give sub-leading contributions and we will not include them explicitly in the following 
expressions. From the group structure there are several possible contractions for matter fields 
ġ with the form of Kahler potential 
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where d, f,d, f denote the corresponding combination coefficients and a is the family index. 
It can be checked that only the first term in Eq. (2.15) contributes in this scenario. Then 
the soft sfermion masses are given by 


FaF Fiste FisFIe 16 
22 16/16 16 16 16* 16 Fis 
emg = ihi- yp +difi Met daft Mt dy fi M2 (2.16) 
after we take into account the normalization factor €? = v?e/M2 of the kinetic term. So we 
have 
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We can see that the typical gaugino mass scale m4/2 is suppressed by a factor e relative to 
the sfermion mass scale. 
The soft SUSY breaking masses for the Higgs potential can be similarly obtained 
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with also contributions from F54 after we add similar terms involving the 54 representation 
Higgs fields. Both contributions are at the same order and the soft SUSY breaking Higgs 
masses are given by MH a ~ (mı 2. Thus we can see that the soft SUSY breaking Higgs 
mass parameters can be at the same order as the gaugino. 

The trilinear terms will also get contributions from both 16 and 16 representation Higgs 
fields. The relevant non-renormalizable superpotential has the form 


W> Yyukawa (T Q Õe)! [C1 (16; & 16i)1010™] 
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Again, here only the first term contributes. From the lowest component VEV and F-term 
VEV of Hı6 and Hyg, we can obtain the trilinear coupling as 


2A, = c12 , vie i6 


E BT Yyukawa - (2.20) 


after we normalize the kinetic term for the 16 representation matter contents. Thus, we see 


that 
A Fs 
Ag = = y E (2.21) 
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is typically of the same scale as sfermion masses. 


The B,-term and p term are given by 


W D (M + ®45)¢10¢10 + kı(M + 45) 
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We will not discuss in detail the doublet-triplet splitting problem in GUT and just use the 
fact that the combination M + (®45) will lead to (by possible tuning or some mechanism) 
proper low energy effective u term. 

The B,-term will be generated after ®54 acquires F-term VEVs 


By = kya / > 54u. (2.23) 


The electroweak symmetry breaking condition can relate different parameters as 


E 4} ond, — tam? pm, (2.24) 
2 tan? 8—1 Ha Hal? 


2B, = sin2B(m%z, + mir, + 2\u|”) . (2.25) 


Although there are many possibilities for the choice of u from electroweak symmetry breaking 
condition, the relation Eq.(2.23) will further constrain the choice of u. Combining Eq.(2.23) 
with the (RGE modified) approximation may, > mi,» we can see that the symmetry breaking 
condition requires |u|? ~ m3, ~ B~ m? J2 for a negative mir. For a positive Mi» it is 
difficult to reconcile a large tan 8 and mir 2 mi, in case of y < mM; Only if the RGE 
running can affect greatly the GUT scale relation may, = mi, can a relatively heavy u be 
possible. The cases of tan 8 ~ 1 which indicate |u|? > [mhl for both signs of mj, can hardly 
be compatible with Eq. (2.23). 

So we see that in our scenario the hierarchy between the GUT scale and the Planck scale 
is used in splitting the SUSY soft spectrum in contrast to the gauge loop factor appeared in 
simply unnatural supersymmetry [14]. 


3 Gauge coupling unification and dark matter constraints 


3.1 Inputs of our model 


In this section we check the gauge coupling unification and dark matter constraints in our 
scenario. From the previous section, we can see the inputs of our scenario at the GUT scale: 


e The gaugino masses with the raito 


Mı : Mə : M3 = —1 : —3: 2 


e The hierarchy between the GUT scale (SO(10) breaking scale) and the Planck scale e. 


e The higgsino mass u which should be at the same sacle as gaugino mass. 


e The parameter tan 8. We scan it in the range 1 ~ 50. 


The sfermion mass mo, determined by mj /2/e, and the trilinear term Ag ~ mo are not 
independent parameters in our model. This scenario predicts different parameter values in 
comparison to the universal gaugino mass scenario [15]. It is well known that the ratios of 
gaugino masses and the corresponding gauge couplings in SUSY are RGE-invariant up to 


d (M; 


If we assume universal gaugino masses Mı = Mə = M3 = My at the GUT scale, from the 


one-loop level 


RGE-invariant ratio 
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and the electroweak scale inputs we can have a mass relation given by 
Mı : My: M3 ~1:2:6 (3.3) 


at electroweak scale. On the other hand, if we assume non-universal gaugino masses at the 
GUT scale, the above gaugino mass ratios at the electroweak scale will be approximately 
changed to 


Mı : Mə : Mz ~ 1:6: —12. (3.4) 


We should note again that the formula Eq. (3.1) can no longer be valid below Mg where heavy 
scalars decouple. So the exact gaugino mass ratio in split SUSY should be obtained by the 
subsequent RGE running from Mg to EW scale. 


3.2 The gauge coupling unification requirement 


We study the gauge coupling unification with the two-loop RGE running of gauge couplings, 
taking into account all weak scale threshold corrections. The relevant analytical results for 
the two-loop beta functions are given in the appendix. It is well known that the two-loop 
RGE running for gauge couplings are scheme independent [23], so we use the MS couplings 
in our study of the gauge coupling unification. Lacking the full knowledge responsible for 
GUT symmetry breaking, especially the mechanism used to solve the doublet-triplet splitting 
problem which is model dependent, we therefore neglect such GUT scale threshold corrections 
in our study. In order to make our calculation reliable, the first step SU(5) GUT scale must be 
much lower than the Planck scale so that the gravitational effects can be neglected. Besides, 
the GUT scale can not be too low, otherwise it will lead to fast proton decay. The constraint 
is [15, 24, 25] 


1.0 x 10'°GeV > Maur > V35acur (6.9 x 10'°) GeV . (3.5) 


In our numerical study we input the central values of gı and g while for g3 we require it in 
the 30 at the electroweak scale. Other inputs at the electroweak scale, for example, the top 
Yukawa coupling y, are extracted from the SM taking into account the threshold corrections. 
Relevant details can be seen in the appendix of [15, 26]. Because of the uncertainty of the 
GUT scale threshold contributions, we adopt the criteria that the gauge coupling unification 
is satisfied when the three couplings differ within the range < 0.005. 
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Figure 1. The scatter plots of the parameter space with gauge coupling unification in case of universal 
gaugino input (upper) and non-universal gaugino input (lower) at the GUT scale. 


Fig.1 shows the result of the parameter space (M2, Ms, p) with successful gauge cou- 
pling unification in case of universal condition Eq.(3.3) and non-universal condition Eq.(3.4), 
respectively. We can see that the gaugino unification gives a stringent constraint on the 
parameter space. From the upper-left panel we can find an upper bound for Mg, which is 
about 10° GeV. Since split SUSY requires Ms > A g We can also obtain an upper bound 
for Mə correspondingly. From the upper-right panel we can find upper limits for u and Ma, 
which are around 100 TeV, independent of the Ms value. However, the constraints for the 
non-universal gaugino scenario are rather mild. From the lower panels, we can see that Mg 
can be as high as 10! GeV while Mz can be 10? GeV. The u parameter, which plays an 
important role in gauge coupling unification, also has an upper bound around 10° GeV for 


both non-universal and universal case. 


3.3 Dark matter constraints 


Now we study the dark matter constraints in our scenario, using the latest dark matter relic 
density data from Planck [27], WMAP[28] and the direct detection limits from XENON100 
[29] and the LUX [30]. The package DarkSUSY [31] is used to scan the parameter space of 
our scenario. In addition to the inputs defined above, we also require 1 < tan < 50. In 
order to calculate the relic density of dark matter, we use the fact that the effects of heavy 
sfermions and heavy Higgs fields almost entirely decouple when Ms = Ma > 5TeV [32]. So 
in our numerical study, we keep the samples which satisfy the GUT constraints and then set 
Ms = Ma = 10 TeV in DarkSUSY to carry out the dark matter related calculations. 

For illustration we set Ms = 30, 50, 100, 200 TeV and show the allowed ranges for 
the parameters Mə, u, tan 8. In our scan, we take into account the current dark matter and 
collider constraints (the details can be found in our previous work [15]): 


(1) We require that the relic density of the neutralino dark matter satisfies the Planck result 
Qpm = 0.1199 + 0.0027 [27] (in combination with the WMAP data [28]). 


(2) The LEP lower bounds on neutralino and charginos, including the invisible decay of 
Z-boson. For LEP experiments, the most stringent constraints come from the chargino 
mass and the invisible Z-boson decay. We require mz+ > 103GeV and the invisible 
decay width [(Z > XoXo) < 1.71 MeV, which is consistent with the 20 precision EW 


measurement [v0" TSM < 2,0 MeV. 


(3) The precision electroweak measurements. We require the oblique parameters [33] to be 
compatible with the LEP/SLD data at 20 confidence level [34]. 


(4) The combined mass range for the Higgs boson: 123GeV < Mp < 127GeV from ATLAS 
and CMS collaborations of LHC [1, 2]. In split SUSY, due to large Mg, log(m=/mé) > 1 
will spoil the convergence of the traditional loop expansion in evaluating the SUSY ef- 
fects of Higgs boson self-energy. So in order to calculate mass of the SM-like Higgs boson, 
we use the RGE improved effective potential(35] which is employed in the NUSSMTools 
package|36] after we set \ = Kk > 0. 


Note that the spin-independent (SI) dark matter-nucleon scattering rate is calculated with 
relevant parameters chosen as [37-39]: roe = 0.023, i = 0:032, = 0.017, N = 0.041 


and ge = S = 0.020. We take into account all the contributions known so far (including 
QCD corrections) in our calculation of the scattering rate. The value of fr, is taken from the 
recent lattice simulation results [40]. 

It is instructive to compare the dark matter constraints for non-universal and universal 
gaugino scenarios. Results for the universal gaugino scenario are taken from our previous 
work [15]. In Figs.2 and 3 we show the samples surviving the constraints (1-4), where the 


green °x’ and red ’A’ denote respectively the samples allowed and excluded by the gauge 
coupling unification requirement. For these results we have the following discussions: 


e We can see from Fig.2 that increasing Mg tends to slightly relax the gauge coupling uni- 
fication constraints in non-universal gaugino scenario, in contrast to universal gaugino 
scenario where increasing Ms tends to spoil the gauge coupling unification. 


e For both scenarios, a strip corresponding to the higgsino dark matter with mass range 
from 1.0 to 1.3 TeV can always survive the combined constraints of gauge coupling 
unification and dark matter direct detection (except Ms > 200TeV in universal gaugino 
scenario which is not preferred by gauge coupling unification requirement). This is the 
well known fact that higgsino at about 1.2 TeV can be a viable dark matter candidate. 


e Outside the strip of higgsino dark matter, almost all the survived points will be covered 
by XEONIT in both scenarios. An interesting exception occurs in non-universal gaugino 
scenario, where a tiny strip at about 50 GeV cannot be covered by XEONIT and such 
a strip enlarges as Mg increases. However, a careful analysis indicates that this strip 
corresponds to a bino dark matter. Although the gaugino mass ratio Mı : Mə : M3 = 
1:6: 12 is no longer valid at the weak scale in split SUSY, the RGE running in general 
will not change significantly the mass ratio for a not too large Mg. So we can estimate 
that this strip corresponds to a gluino below 700 GeV. The current preliminary limits 
on gluino mass using 20 fb! of 8 TeV data are Mj = 1350 GeV ([43]) and Mj = 1200 
GeV ({44]) assuming a massless neutralino for mini-split SUSY [45]. So this tiny strip 
should have been ruled out by the LHC data. 


So we can see that a neutralino dark matter below about 1.0 TeV will be fully covered 
by XENONIT for both non-universal and universal gaugino scenarios in split SUSY. We 
checked that our numerical calculation results are not sensitive to the sign of u, except that 
for the minus sign a very small part of parameter space (called blind spots) can survive all 
the constraints including the future XENONIT limits [42]. 


4 Conclusion 


In this work we proposed to achieve the hierarchy of sparticle spectrum in split SUSY from 
the gap between the GUT scale and the Planck scale. We built a split SUSY model (which 
gives non-universal gaugino masses) with proper high dimensional operators in the framework 
of SO(10) GUT. Based on a calculation of two-loop beta functions for gauge couplings (taking 
into account all weak scale threshold corrections), we checked the gauge coupling unification 
and dark matter constraints. We found that our scenario can achieve the gauge coupling 
unification and satisfy the dark matter constraints in some part of parameter space. We also 
examined the sensitivity of the future XENONIT experiment and found that the currently 
allowed parameter space in our scenario can be covered for a neutralino dark matter below 
about 1.0 TeV. 
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Figure 2. The scatter plots of the (u, M2) parameter space satisfying constraints (1-4) including 
dark matter relic density. The green ’x’ (red ’A’) can (cannot) achieve the gauge coupling unification. 
The left panel is for the non-universal gaugino scenario proposed in this work while the right panel is 
for the universal gaugino scenario studied in our previous work [15]. 
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Figure 3. Same as Fig.2, but showing the spin-independent cross section of dark matter scattering 
off the nucleon. The curves denote the limits from XENON100(2012) and LUX as well as the future 
XENONIT sensitivity. 
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A The two-loop beta function for gauge couplings 


In this work, we adopt the method in [46, 47] to calculate the two-loop beta functions for three 
gauge couplings in the split SUSY, taking into account the weak scale threshold corrections. 
Our result is in agreement with [11, 13]. The notations for the two-loop RGE can be found 
in the appendix of our previous work [15]. The one-loop beta functions for the three gauge 


couplings are given by 


bs = —76( — Mz) +28 — M3) + 20(u — Ms) , (A.1) 
by = alu Mz) 4 “Au Myy) + Au — My) + oqu —Ms); (A2) 
bi = Hofu- Mz) + 20(u— Ma) + 70l- Ms), (A.3) 


with the step function defined as 


1 x > 0; 
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The 2-loop RGE for SU (3)e, SU (2)z, U (l)y gauge couplings (g3, g2, g1, respectively) are 
given by 
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with the U(1)y normalization g? = 3(gy)?. 

The two-loop beta functions for gauge couplings are given by 
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20 4? 20’ 4 
B = (2, * 0)6(u— max(Mi,Mg)) + (u— Ms) -Ż,-=,0. (AL) 
20’ 4’ me 20° 4’ 


The one-loop renormalization group equations for Yukawa couplings below the Mg scale can 


be written as 


d 3 
16m? = he =h" -eto + fT + 8,51 Hes (arint = niin) ,  (A.12) 
16r? p =h’ [-3eto? + chT + ch S1 + ch, S2 + : (a = nin) ,  (A.13) 
2 d e e e 2 e e e 3 etpe 
167 a = h? |—3c g; + epTT + cg, S1 + Ch, S2 + 5h hfj, (A.14) 
with 
1 3 
T = Tr(3hth" + 3AA? + Ath’), Si = 5 [(Gu)” + Ga)"] , S2= 5 (Ga + gå) - (A-185) 


Above Ms, we will recover the MSSM result and the one-loop RGE for Yukawa-type inter- 


actions in the superpotential are well known to be 


161?" =)" |—2ci'9? + 3Tr (AM A") + BAMTAY + adtad] ; (A.16) 
16? 4 = \4 | 2046? fT PCBNT AS + ATAS) + AYTAY + aatar] ; (A.17) 
d 
16m? = A° =)? |-2¢F9? + Tr (BAT At + ATAS) + Bastat] ; (A.18) 
with 
13 3 8, 4 738 9 3 
Ui (le Ss s es Pe E = Ce . A.19 
Ci (aga? a)? Ci (357 93) Ci (597 99) ( ) 


The one-loop Yukawa couplings for c¥, ef, c; are calculated to be 


cl 17 3 8 2 39 
i 60 4 3 20 4 
cf | =O(u—Mz)| $ #3] +0u-—Ms)| $30 (A.20) 
cf 7 70 3 70 
All terms are set to zero above Mg. Besides, we have 
Ch Cs, CS, 100 001 
ch ch ch, | = O(u— Mz) | 100 | +0(u-max{ M, Mg) | 001] , 
Cr CS, Co, 100 001 
010 
+ 0(u — max(M,,M;,)) | 010 (A.21) 
010 


— 13- 


and all set to zero above Mgs. 


One loop RGE for gaugino couplings ĝu, Ji, Ja, J4 below Msg are given as (with the gaugino 


relation Mı < M2) 
e Between [max(M2, M;,),Ms], the RGE for u,a are given by 


d 7 ee ae ee ee ere 
16n? 7 Gu = —3Guct g; + qiu — ziela + gõuu + GaGa + Gu(T + cs, $1 + es, 2) 
a - EO ee ee ee eee ee 
167? —Gq = —3gaclg? + 203 — =GaGn + gag + Juuda + Ga(T + cs, $1 + 5,92) 
dt 4 2 4 
with the coefficient 
T wee 
G = (= —,0), Cs, = CS2 = 1. 


Below max( M2, Mẹ), the coupling are switched off. 


e Between [max(M2,M;,),Ms], the RGE for gq are given by 


d. ye Bole a Op har. Bhs See ee 
16r? Gu = —8Guet oF + Gie + giuda + gludu + 3GGaGu + G(T + es, $1 + es,$2) 
1672 iz _ 39,0492 4 3213 | 32) 72 | 3 1 a2 4 gal a õa + 5L(T Si +¢¢,8 
q Jači Gi + 79d + gau + qala + 39uIuGa JalT + eg, S1 + cs3 82) 
with the coefficient 
ud 7% 3 O O 
i 7 (55° po Cs, = Cs, = 1. 


Between [max( M, Mẹ), max( M2, M;;)|, the RGE reads 


d. n Jea er es 
16r? Ie = BHU + He + FLIE + GT + 05, 81 + ¢5482), 
a te a ae 
16r’ Gi = —39)0%97 + zid + 59419 + gh(T + cs, S1 + cs3 S2), 
with the coefficient 
~u,d 3 
a= (o rO cg, = 1, cs, =0 


Below max( Mı, M;;), the coupling are switched off. 
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